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Abstract 

The most commonly employed formalisms of radiative energy loss have been derived in the 
high- energy approximation. In its present form, it is reliable only for the medium modifications 
of inclusive particle spectra. Modifications to this formalism are expected to be important for less 
inclusive measurements. This is especially relevant for reconstructed jets in heavy-ion collisions, 
which are becoming available only recently. We present some ideas to overcome this limitation. 
Specifically, we show an implementation of radiative energy loss within a jet parton shower. This 
implementation has been done within the PYTHIA Monte Carlo event generator. We present the 
publicly available routine Q-PYTHIA and discuss some of the obtained physics results. 



1. Introduction 

Reconstructed jets, a long-awaited probe of medium properties, are starting to become a 
reality in heavy-ion collisions |1, 2]. The expectations at the LHC are very high for these new 
probes jsl]. On the theory side, the description of a reconstructed jet as a probe of the medium 
has already some history. However, analytical approaches are very limited to deal with such 
a complicated signals — this is nothing special in heavy-ion collisions but common to other, 
more simple, systems as well. In this conditions, Monte Carlo codes to describe the jet structure 
with the needed degree of exclusivity become essential tools for both the phenomenological 
and the experimental analyses — where corrections based on these codes are frequent. Several 
implementations have been presented in the last year, each of them trying to address a particular 
question of the whole problem. 

On the theoretical level, jets are hard processes which can be described in perturbative QCD. 
The factorization theorems allow to separate long- from short-distance parts in the cross section. 
The long distance parts encode all the non-perturbative contributions to the process. An inte- 
gral part of the cross sections is how these long distance parts evolve with virtuality, which can 
be computed using perturbative techniques by resumming the relevant splitting probabilities en- 
hanced by phase space factors. The so-called final-state radiation leads to a description of the jet 
structures. The first and most naive expectation in the case that a medium is formed in heavy-ion 
collisions is that this radiation pattern is modified. Indeed, most of the successful phenomenol- 
ogy of high-/77 particles at RHIC relied on the radiative energy loss as the main source of the 
effect — see e.g. Ref. for a recent study of experimental data. What we present in here is a 
particular implementation on how this medium-induced radiation can be included in a final-state 
parton shower and how the corresponding jet structures are modified. 
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2. The modified splitting probability 



The gluon radiation spectrum in the presence of a medium, when computed including the 
important interference with the vacuum contribution, can be split into a medium-induced contri- 
bution plus a vacuum contribution [SJ 



dojdkj dojdkj dcodkj 

This factorization has been shown at the level of the one-gluon inclusive distributions. Here, 
however, we are interested in the exclusive one-, two-, three- etc gluon distributions. In the vac- 
uum, the different interference and virtual terms are known to a given accuracy and the dominant 
contribution for a multi-parton emission is identified to be given by the ordered (in kr or angle) 
emission of subsequent gluons, the non-ordered contributions being subleading. A similar com- 
putation has not yet been performed for the case of the medium, where, in fact, one would expect 
that other ordering variables related with the space-time structure of the shower could become 
relevant — this space-time structure is irrelevant in the vacuum where no length scale external to 
the shower is present. 

In our implementation, we use a mixed approach, in which the splittings are generated as in 
the vacuum — with a similar ordering variable — but corrections for the space-time evolution are 
introduced at every splitting. This allows us, in particular, to treat medium- and vacuum-splittings 
on the same footing (0, 0] . In particular, we define a medium-modified splitting probability as 

Ptotiz) = Pvac(z) ^ PUZ) = Pvac(z) + AP(z, f, q, L, E). (2) 

By matching with the vacuum case, the corresponding modification to the splitting probability is 
written as 

2nk^, dl"""^ 

AP(z, t, q, L, E) = i J (3) 
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where dl"^'^'^ I dcodkj is given by the medium-induced gluon radiation used previously for RHIC 
phenomenology [5]. This spectrum depends on two parameters, the medium length L and the 
transport coefficient q. The latter of these quantities encodes all possible information about the 
medium properties as temperature, density, etc. in a single parameter, with the interpretation of 
the average transverse momentum squared that the gluon gets per mean free path in the medium. 



3. Implementation and results 

The medium-modified splitting probabilities (|2] O are implemented in a modification of the 
standard PYTHIA [8] routine for final-state parton branching PYSHOW. There, the vacuum split- 
ting functions are supplemented with the medium term — see |6] for details. The modified 
routines along with auxiliary routines to compute the medium contribution to the splitting prob- 
ability and geometry are released for public use and available at the Site |9]. As a short name, 
we call this implementation Q-PYTHIA [6]. 

The additional medium-induced splitting probability is expected to produce: 

1 . Softening of the spectra, leading to an apparent energy loss of the leading particle in the 
jet 
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Figure 1: (Left) Jet shapes in pp collisions at ^JSNN = 5.5 TeV using tlie anti-kr recombination algorithm for different 
quenching parameters. (Right) Ratio of medium to vacuum fragmentation functions for all hadrons, for different gluon 
energies Ejg, = 100 GeV. Figures from Ref. 

2. Increase of the intra-jet multiplicity 

3. Broadening of the jet angular structures due to the larger typical angle of radiation in the 
medium when compared to the vacuum. 

These features are observed in the physics results from our implementation Q-PYTHIA l'6'l. In 
Fig. [T]the angular energy broadening and the softening of the fragmentation function is plotted 
for different geometries and values of the transport coefficient. It is worth noticing that a very 
large energy broadening could induce a bias in the jet reconstruction while (due to the softening 
of the spectrum) the study of intra-jet particle angular distributions is a good observable for 
medium effects ifiol . 

4. Usage of the routines 

Q-PYTHIA is standard PYTHIA with a modified final-state parton shower routine. So, all 
the features of this Monte Carlo are unchanged except for the evolution of the jet after the hard 
scattering — we refer to the standard manual ^ for usage. 

The medium-modifications are introduced in Q-PYTHIA through two geometry routine^ 
QPYGIN to define the initial position of the hard scattering and QPYGED which returns the medium 
properties that a parton traversing the medium will measure. In order to allow for a flexible im- 
plementation, the geometry of the medium is not fixed but can be defined by the user — examples 
of simple geometries are provided and more of them will be available in the future. The medium 
enters in the splitting probabilities through only two variables, the transport coefficient, q, and the 
in-medium path length, L. In fact, the actual variables used in the implementation are a combi- 
nation of these two: qL and = qL^/2. For a generic medium in which the transport coefficient 
is defined locally q(x, y, z, t) one usually defines these quantities as the effective integrals 

cj'/^(x,y,z,T,f5„p,,f5,)^ f d^m). 



Geometry here is to be understood in a generic way as the space-time extension of a medium, including expansion, 
etc. In fact, what is needed is a transport coefficient defined locally at every point of the space and every time after the 
collision. 
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{qLrHx,y,z,r,p,,liy,p,)^ jd^qi^ (4) 

where x, y, z and r are the initial position and time to be considered (the initial production point 
or the point of last splitting) and yS,, / - x,y,z, the components of the corresponding trajectory 
three-velocity in units of c. The integrations in Eq. (|4]l are to be performed along the linear 
trajectory defined by these inputs and parametrized by ^. 

Once the initial point of the jet evolution is fixed by qPYGIN(XO,YO,ZO,TO), Q-PYTHIA 
will start the jet evolution by computing the probability of splitting and assigning it a space-time 
position computed by the formation time of the secondary parton, estimated as 

2(aj 

^form ~ "^2~ 

This time defines then the subsequent position (x, y, z, t) where the next splitting will start — Q- 
PYTHIA calls then QPYGED (X , Y , Z , T , BX , BY , BZ , QHL , 0C)| which returns the values of [qL]"^^ 
and a//-^ for this particular path. In this way, the formation time is introduced at every splitting 
allowing to account for the space-time evolution of the jet structure. So, the standard way of 
implementing a geometry defined in the routines QPYGIN and QPYGEO simplifies to: i) defining 
a local transport coefficient q{x, y, z, t) — for example by assuming q oc e^/'^ and e the energy 
density field computed in a hydrodynamical simulation; ii) computing the linear integrals (|4]i for 
each path defined by the inputs (x, y, z, t) and (/3_^,/3y,P;). 
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^When the splitting is the first one, (x,y,z,T) will coiTespond to the initial values (j:o,>'o,zo,^o)- 
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